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The nucleon-nucleon potential is analysed using the 1/Nc expansion of QCD. The NN 
^ , potential is shown to have an expansion in and the strengths of the leading order 

2 . central, spin-orbit, tensor, and quadratic spin-orbit forces (including isospin dependence) 
are determined. Comparison with a successful phenomenological potential (Nijmegen) 
^ , shows that the large- A'^c analysis explains many of the qualitative features observed in the 
■ nucleon-nucleon interaction. The l/N^ expansion implies an effective Wigner supermulti- 
plet symmetry for light nuclei. Results for baryons containing strange quarks are presented 
in an appendix. 
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1. Introduction 



The two-body nucleon-nucleon interaction is the basic ingredient that is used in cal- 
culating the properties of nuclei. There are a number of phenomenologically successful 
models for the interaction, typically constructed using meson exchange contributions. Un- 
fortunately there is no direct connection between these models and the underlying theory 
of the strong interactions, QCD. Until recently, the only way that QCD and the physics 
of nucleon interactions could be rigorously related has been through symmetry arguments. 
By making use of symmetry and effective field theory, one can calculate the low energy 
dynamics using a small number of parameters that are fitted to the experimental data. 
Such theories can be quite predictive, but still remain somewhat remote from QCD — one 
would like a priori arguments for the sizes of the phenomenological parameters. 

Recently, there has been significant progress in better understanding the implications 
of QCD for hadronic physics by exploiting the "hidden" expansion parameter of QCD — 
1/Nc, where A'^c = 3 is the number of colors In the large Nc limit, one finds that 

meson-baryon interactions respect an SU{4) spin-flavor symmetry — the same symmetry 
found in the nonrelativistic quark model quark model It has also been shown that 

the 1/Nc expansion can provide information about the nucleon-nucleon potential 0. In 
particular, ref. analyzed the central potential for A^A^ scattering and showed that the 
expansion gave a qualitative understanding of its spin and isospin structure. 

In this paper we extend the analysis of to include the entire A^A^ potential. Naively, 
one might think that the Nc ^ oo limit is not relevant for analyzing nuclear physics. 
Nuclear matter forms a classical crystal at Nc = oo, and so there must be a phase transition 
between Nc = 3 and Nc = oo. While the 1/Nc expansion is not reliable for studying bulk 
properties of nuclear matter, it does allow one to analyze the spin and isospin dependence 
of the nuclear force. One expects that the symmetry properties of the A^A^ interaction will 
be independent of the phase of the many body groundstate. 

The general form of the potential for elastic, nonrelativistic A^A^ scattering is 

Vnn = + y>i ■ (72 + V^s^ ■ S + VSSi2 + V^Qi2 

(1.1) 

+ (Fo' + ■ ^2 + Vlsh ■ S + V^Si2 + ViQi2) ri ■ r2 , 

where 

5*12 = 3(Ti • f (72 ■ f - (Ti ■ 0-2 

(1.2) 

Qi2 = i{(^i-L),(a2-L)} . 
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TABLE 1 



Isospin 


Vo 




Vls 


Ft 




1 1 

ri ■ T2 


1/iVc 


1/iVc 


1/iVc 
1/iVc 


1/iVc 


l/iV3 
1/iVc 



The four terms Vq, constitute the central potential, while V^, Vlg and Vq are the tensor 
interaction, the spin-orbit interaction, and the quadratic spin-orbit interaction respectively; 
the ten functions V^ can in general be velocity dependent. The main result of this paper 
is that the strength of the ten functions V^ can be determined in the 1/Nc expansion, and 
are as given in Table 1. As is apparent from this Table, the actual expansion parameter 
is not 1/Nc but Thus even though the actual value Nc = 3 is not very large, an 

expansion in l/N^ can be quite predictive. 

The organization of this paper is as follows: in §2 we briefly review general properties of 
baryons in the large Nc limit. In §3 we derive the results given in Table 1. These results are 
compared in §4 with the "Nijmegen potential" of references — a phenomenologically 
successful model of the A^A^ interaction; we show that the the hierarchy of Table 1 is 
evident in A^A^ phenomenology. §5 extends the discussion of ref. concerning how the 
Wigner supermultiplet symmetry might arise in light nuclei as a consequence of the 1/Nc 
expansion. This is followed by conclusions, and an appendix in which we extend our 
analysis to hyperon interactions (z.e, including the s quark). 



2. The Large A^^ QCD Analysis 

The large Nc limit is deflned by taking the number of colors Nc of QCD to be large 
while simultaneously rescaling the QCD coupling as (7 — > g/ ^/Wl, keeping Aqcd flxed |T^ . 
The 1/A^c expansion has proven to be a powerful tool for analyzing baryon properties, 
since baryon structure simplifles considerably in this limit. Antiquarks in the baryon 
are suppressed, and as baryons consist of Nc quarks interacting with 1/A^c strength, the 



Hartree approximation becomes exact in the large Nc limit |]Tl| . Although one cannot solve 
the Hartree equations due to the nonlinearity of glue interactions, one can nevertheless 
determine a number of useful properties of the spin and flavor properties of the baryons 
and their interactions. 
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Fig. 1. A QCD contribution to H leading in 1/Nc. This diagram can be described 



in spin-flavor space as a product of three of the 1-quark operators given in eq. ( ^.1\) . 



To analyze the flavor and spin structure of the Hartree Hamiltonian in the case of two 
hght flavors, it is convenient to use as an operator basis the one-quark operators of the 
quark model 



^^ = 9^y9, /" = 9^y9, G^'' = q^^q, (2.1) 

where q = (w, d) and are the creation and annihilation operators for the u and d 
quark flavors, and a*, are the standard SU{2) Pauli matrices acting on spin and isospin 
respectively. The q and operators do not carry color, and are bosonic. The Hartree 
Hamiltonian can then be constructed as monomials of these operators. An important 
result from large Nc QCD is that the Hartree Hamiltonian takes the form 



n s,t 



where the operators {S, /, G} are given in eq. (|2.1j), the coefficients v are 0{1) functions 
of momenta, and we have suppressed isospin, spin, and vector indices which are contracted 
such that H is rotation and isospin invariant. An example of a contribution to H is pictured 
in fig. 1. It is important that although we make use of the quark model operator basis, 
eq. ( p.2|) makes no assumption about the validity of the quark model; the quark model 
operators are a representation of the spin-fiavor Clebsch-Gordon coefficients, and provide 
an efficient way of doing group theory computations. A Skyrme model basis, for example, 
would have worked just as well . 



The lowest lying eigenstates of H with baryon number B = 1 have isospin / and spin S 
satisfying / = 5'=^,|,|.... The first two states can be identified with the and A; the 
higher spin states do not exist for Nc = 3. To leading order, this baryon tower is degenerate 
with mass M ~ Nc- One can show that it transforms as an irreducible representation (the 
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totally symmetric A^c-index tensor) of an approximate SU{4) spin- flavor symmetry; this 
is the symmetry under which the quark operators m|, w|, (if, (i| transform as the 
four-dimensional fundamental representation |jl|,|13[. For Nc = 3 the {A^, A} spin states 
transform as the 20 dimensional representation of SU{4), familiar from the quark model. 

Additional information can be obtained by considering matrix elements of operators 
between baryon states B and B' restricted to have S — I ^ 1. For example, matrix 
elements of the basis operators (^3) satisfy 



{B'\ S/N, \B) ~ {B'\ i/N, \B) ~ , {B'\ G/N, |S) ~ 1 . (2.3) 

Matrix elements of many body operators can be analyzed as well, using various relations 
among powers of the basic operators S\ and This allows one to greatly 

reduce the number of linearly independent terms in ( p.2|) at a given order in 1/Nc- Using 
these techniques, one is able to show that the matrix element of a general n-quark operator 
O^Pl with S = 0, isospin / and spin S, is of size (see also for a derivation) 

{B'\ dg/iV- \B) < l/N\J-^\ . (2.4) 

The fact that the operators with the largest matrix elements have I = S was first observed 
in the Skyrme model [0, and is known as the "/t = Jt rule". 



Eqs. ( p.2|) -( |274D are the central results behind the large A^c analysis of baryons. One 
consequence is that the mass splittings in the baryon tower (e.g, between N and A) 
are of size P JTl]] . This result is in good agreement with the real world where the 

ratio R = (Ma - Mn)/{Ma + Mn) = 0.13, while the large prediction at A^c = 3 is 
R ~ 1/A^c — 0.11. Consequences of eq. ( |2.3|) for the A^A^ interaction are explored in the 
next section. 



3. The Nucleon-Nucleon Interaction 

There are two independent three-momenta for baryon-baryon scattering in the center 
of mass frame, which can conveniently taken to be 

(l=Pin -Pout , ^ =Pm +Pout " (3-1) 

These momenta are to be considered independent of Nc in the 1/Nc expansion. To leading 
order in 1/A^c, the entire baryon tower is degenerate and = \Pout \ ^or elastic scattering. 
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up to 0{1/N^) corrections, and so g ■ /c = to the same order. The general baryon-baryon 
interaction potential is then a matrix 

V{q, k) = (p„^i, 7; -p^^t, 5\ H \p^^, a; (3) (3.2) 

where H is the Hartree Hamiltonian (|2.2| ) and a, ... ,5 denote internal quantum numbers 
of the baryons, such as spin, flavor and particle type {e.g, N or A). Throughout this paper 
we will deflne the A^A^ potential as the above matrix element restricted to the space of 
nucleons. We do not consider second order effects due to virtual A's, etc. 

There are two ways that 1/Nc factors can suppress terms in the potential. The flrst 
arises from spin-flavor structure and the powers of 1/Nc in eq. ( |2.3| ). The second source 
of suppression arises in velocity dependent interactions arising as relativistic corrections. 
Since the nucleon velocity equals p/M ~ l/A^c, each power of velocity is equivalent to a 
1/Nc suppression. In the nonrelativistic limit for baryons, a t-channel meson exchange 
contribution to V is only a function of q. A w-channel contribution is only a function of 
k, and can be expressed as an exchange potential. Relativistic corrections allow a single 
meson exchange contribution to y to be a function of both q and k. Meson exchange in 
the t-channel is then a function of q and k/M, with each power of k being accompanied by 
one factor of M. Similarly, w-channel meson exchange is a function of k and q/M. This 
shows that if a general velocity dependent potential is expanded in a Taylor series in k and 
q, a term of the form q^k^ is suppressed by 

l/NJ^ , n = Min(r, s) . (3.3) 

Combining this source of 1/Nc suppression with eq. (|2.3|) will allow us to determine the 
size and spin-flavor structure of the dominant terms in the potential V. 

An A^A^ interaction at the QCD level gets contributions from complicated processes, 
such as pictured in flg. 2. Each of these contributions can be expressed as a tensor function 
v{q,k) contracted with 1-quark operators 5"*, and G*" which act on either of the two 
nucleon states. The coeflicient function Vstn in eq. (|2.2|) and the operators S'^ and G*" must 
combine to be invariant under rotations. Our analysis is simplifled by flrst expanding Vstn 
(and hence V) in multipole moments; the following subsections are organized accordingly 
as AL = (the central force), AL = 1 (spin-orbit force), and AL = 2 (the tensor and 
quadratic spin-orbit forces). 
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Fig. 2. An example of a contribution to the NN interaction from the 3-quark 
operator pictured in fig. 1. This diagram can he described in spin-flavor space as 
a single 1-quark operator acting on the first baryon Ni, and two 1-quark operators 
acting on the N2 line. 



3.1. AL = 0; The Central Potential 

This is the case analyzed in ref. 0]. The central force can be written as a sum of 



products of 1-quark operators as in eq. (|2.2|) , where the operators act on either the A'^i or 
N2 nucleon states, and the coefficients Vstn are general scalar functions of |g| and |fc|. It 
follows from eq. (|2.3D that the leading contribution will have no S'^/Nc or I°'/Nc operators, 
since each of these implies a 1/Nc suppression; instead it will consist solely of powers of 
G'^"'/Nc. By rotational symmetry, since the v coefficients are scalars, the G*" operators 
must be contracted to form spin invariants. Similarly, isospin symmetry implies that the 
Qia j^^g^ contracted to form isospin invariants. From these constraints, it is possible to 
show that at leading order in 1/Nc, the most general form for the central potential is ^ 




^central =N,}Vn\ ' , , (3.4) 



where Gi ■ G2 = G^^T- general, the coefficients Vn are functions of both |A;|^, 
and obey the rule eq. (pi. 3]). One can further restrict the powers of Gi ■ G2 in eq. ( |3.4[ ) to 
be completely symmetric in the Gi indices, and in the G2 indices, before the two sets of 
indices are contracted. 

It is straightforward to verify that ( |3.4| ) is the most general form of the leading order 
AL = potential. We have argued that it can only involve powers of the G*" operators, 
on the basis of eq. ( |2.3| ); what must be shown is that the indices are contracted as above 
in eq. ( |3.4| ). By the operator reduction rule ^ any terms in which two indices of &'^G^^ 
(where both G's act on the same baryon) are contracted with each other by S^^ , S""^, e^^^ 
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or e"^'^ can be eliminated in favor of terms with fewer powers of G Thus the only allowed 
invariants are obtained by contracting the indices of G^" with those of Gj", as in ( |3.4|) . 
More complicated contractions, such as 



Gf Gf (3.5) 



can be written as 



Gi -GajVcf Gf 



Gi ■ G2) ^ - e"^' Gf Gf e^'^" Gl" Gi 



(3.6) 



The term with two e symbols can be reduced to Gi ■ G2 using the relations in 0, so 
that all contractions of Gi with G2 can be written as powers of Gi ■ G2- One can also 
restrict the indices on powers of Gi and G2 to be completely symmetrized, since terms 
antisymmetric in the indices can be eliminated using the operator identities. The series in 
Gi ■ G2 terminates after Nn terms, because an operator with more than Nr quark fields 
a^ting „. a single baryo. ca. be ree.pre.sed tel of operate, with < ,L Helds.B 
There are no 1/Nc corrections to ( |3.4| ), except through 1/Nc dependence in the un- 
known coefficients Vn- AH terms in the 1/Nc correction to ( |3.4| ) are arbitrary polynomials 
in Gf2 5 with one factor of SI 2 or /"2- It is easy to check that all such terms have the 
wrong time-reversal properties to contribute to the baryon-baryon potential. Thus the 
first correction to (|3.4|) contains a factor of Si ■ S2 or Ii ■ I2 and is of order 



Eqs. ( |2.2D , ( |2.4| ), and (|1.1| ) define the A^c-counting for the central potential. Large 
Nc QCD implies that the central potential is of order A''^, but is determined by only two 
independent functions instead of four (at leading order in 1/Nc): 

Vo^{r)r^Nc, V}{r)r^Nc, (3.7) 

while 

V^ir)r^l/Nc VoHr)-l/Nc. (3.8) 

As was noted in and will be discussed in §5, the above relation implies that the central 
potential obeys an effective Wigner supermultiplet symmetry. 



^ An easy way to see this is to normal order the operators. A normal ordered product involving 
with more than Nc quark operators on a baryon vanishes, which gives the desired identity. 



7 



3.2. AL = 1: The Spin- Orbit Potential 

The AL = 1 baryon interaction amplitude contains the spin-orbit coupHng term; 
it is obtained from the general Hartree Hamiltonian eq. (|2.2|) by restricting attention to 
terms for which the coefficient v transforms as a vector under rotations. It follows that 
the 1-quark operators multiplying the v coefficients must be combined to transform as a 
(1,0) representation under SU{2)spin x SU{2)isospin- From eq. ( PT^ ) we have seen that 
to contribute at leading order in 1/Nc, an n-quark operator must be a polynomial in the 
G's alone. However, one cannot make a (1,0) operator with the correct parity and time 
reversal properties purely out of G"s. The spin-orbit force is suppressed relative to the 
central force, and is an arbitrary polynomial in G's, with one factor of S or /. The general 
form of the AL = 1 amplitude is 




Vls 



n=0 



Ltj^ (j"2 i~>\ "r 'J'2 "J"! i I 1 ' 2 



This can be derived by arguments similar to those in the previous subsection. Time reversal 
and parity invariance requires the coefficients „ in ( ^^) to be proportional to (g x fc) 
times an arbitrary function of and /c^. In position space, a contribution of the form 
U (g^) (g X fc) ■ (5^1 +6^2) is of the form (Vf/ (r) x k) ■ S, which is the usual spin-orbit force. 
There is a hidden suppression factor of 1/Nc (which follows from eq. ( |3.3| )) in the spin-orbit 
force which is not manifest in eq. ( |3.9| ), since the AL = 1 interaction necessarily involves 
both q and k. 

The Wigner-Eckart theorem implies that there are only two distinct operators when 
the expression ( |3.9|) for the AL = 1 amplitude is restricted to the nucleon sector. These 
are the two spin-orbit terms appearing in eq. (pTl]). Thus we find 

VEsir) ~ 1/iV, , Vlsir) ~ 1/iV, . (3.10) 

The spin-orbit force is 0{1/N^) in strength relative to the central force, and it is of 
comparable strength in the two isospin channels. 
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3.3. AL = 2: The Tensor and Quadratic Spin- Orbit Potentials 

The AL = 2 amplitude is obtained by requiring that the coefficients v in eq. 
transform under rotations as AL = 2. The leading order amplitude is a polynomial in the 
(j's that transforms as = 2, / = 0. One can obtain an amplitude that does not violate 
the If = Jt rule on each baryon line by combining / = 5" = 1 amplitudes on each baryon 
into total / = and total AL = 2. The general form of the leading order amplitude is 

n=0 ^ \ ^ / 

where the coefficient f is a symmetric traceless tensor that depends on g and k. Time 
reversal invariance requires the coefficients to have the form 



where is a scalar function of g^, fc^ and (g ■ /c)^. 

If one restricts the interaction eq. (|3.11|) to the nucleon sector, one gets 



= A^c Vn Ti ■ T2 (^q ■ (Ti q ■ a2 - ^q^ ai ■ 



(3.12) 



Terms with n > 1 in eq. ( 3.11 ) can be dropped, because two spin-1/2 nucleons can only 



give non-zero matrix elements for operators with spin < 1. Comparing with eq. ( |1 . 1| ) , we 
see that 

~ N^. (3.13) 

The other term in the tensor potential, V^^, has |/ — 5*1 = 1 at each nucleon line, and so 
by eq. (U) 

Vt ~ 1/Nc. (3.14) 
A similar and straightforward analysis for Vq gives the results listed in Table 1. 

3.4. The NN potential and the A 

One flaw in our discussion that should be eventually improved upon is the treatment 
of the A. The Hartree Hamiltonian (|2.2| ) implicitly acts on the entire I = S baryon tower, 
including both nucleons and A's, all of which are degenerate in the Nc —>■ 00 limit. In our 
discussion of the NN potential, we have simply projected H to the nucleon sector. A more 
sophisticated treatment would be to integrate the A's out of the theory (keeping track of 
the 1/A^c mass splitting) and to construct an effective theory for nucleons alone. This is a 
subtle analysis (see, for example, fl^) and beyond the scope of this paper. 
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4. Comparison of large Nc QCD with a phenomenologically successful model 

Our large- A^c results for the general nucleon-nucleon potential of eq. are displayed 
in Table 1. For two flavors we have found that the strongest A^A^ interactions are the central 
force terms Vq and VJ , as well as the tensor force Vf, all three of which are ~ Nc- The 
remaining contributions to the A^A^ potential, with the exception of Vq, are relatively 
suppressed by ~ 0{1/N^). Finally, the isospin invariant quadratic spin-orbit force Vq 
is suppressed by ~ 0{1/N^) compared to the central potential, as it is both an / 7^ 5" 
interaction, as well as being a second order relativistic effect suppressed by The 
results we have derived are consistent with the It = Jt rule, but are more general. They 
are true in QCD in the 1/A^c expansion, and make no assumptions about the origin of the 
A^A^ interaction as being, for example, due to one meson exchange. 

The results can be directly compared with nuclear potential models in momentum 
space. A particularly simple phenomenological model to compare with is the meson 
exchange model "Nijmegen potential" of references 1^,0. In this model, the A^A^ po- 
tential is approximated in momentum space by a sum of Yukawa and Gaussian interac- 
tions times powers of momenta divided by masses, contracted with the spin and isospin 
Pauli matrices. The Yukawa potentials correspond to one-particle exchange of both real 
mesons(7r, 77, 77', p, w, ^, ao, /o) and an "effective meson" (e), while the Gaussian potentials 
are labelled by P, /2, /2 and a2- The motivation for this form of the potential is unim- 
portant here; it provides a phenomenologically successful parametrization for the A^A^ 
potential that can be compared with Table 1. The Nc dependence should appear in the 
relative strengths of the potentials, and the 1/M factors that appear when the potential is 



decomposed as in eq. (IT). The strength of the contributions to the Nijmegen potential 
are simple to evaluate, since they are presented explicitly in momentum space, and we can 
treat all momenta and meson masses as ~ 1 in the 1/A^c expansion. The Nc dependence 
must then reside in the strengths of the couplings used in the Nijmegen potential, as well 
as the explicit factors of the nucleon mass that appear in the formulas of ref § . One flnds 
for the strength of the various terms in the potential 

yi yi 9jo_ 9io9ii 9n_ 

" ^ ^ ^ M2' AM ' A2' ,^ ^, 

(4.1) 

yi 9jo_ 910911 
M2 ' AM ' 

yi 9jo_ 9I09I1 9'ii 
~ M4' AM3 ' A2M2 ' 
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where / = 0, 1 correspond to the 1 ■ 1 and ti ■ T2 isospin structures respectively, M is the 
nucleon mass, and A is a strong interaction scale characterizing the derivative expansion 
(denoted A4 in ||]). The parameters gis correspond to the coupling constants of the 
model with t-channel (isospin, spin)= (/, 5") in the nonrelativistic limit; in particular, 
the scalar coupling gs and vector couplings gy and fy of ref. [§] are given by gjo, gjo 
and gii respectively, where / is the meson isospin. (The pseudoscalar contributions are 
parametrized differently in ^ and are mentioned below). As far as the Nc scaling goes, 
M ~ Nc, while the A and the masses of the exchange mesons are all ~ 1. In eq. ( [4.1J ) 
we have omitted dimensionful quantities that do not scale with such as the meson 
propagators l/{q'^ + m?). By comparing the expressions in eq. ( f4.1| ) with our results in 
Table 1, one sees that they are consistent provided that the couplings gis scale with A^^ as 



9is oc 



This Nc scaling can be compared with the numerical values given in ref. [^. In fig. 3 we 
have plotted the couplings determined numerically in ref. f^, rescaled by their value for 



fp. Since fp is a, gn coupling, eq. ( [4.2|) implies that the leading large- A'c prediction for 
the ratio is 

g,s f 1, if|/-S| = 

= if|/-si = i ■ 

As can be seen from fig. 3, there is good qualitative agreement between the large- A'c 
prediction (|4.3| ), and the gis values used in the Nijmegen potential. Omitted from fig. 3 are 
the pseudoscalar couplings, which are not readily compared with heavy meson couplings, 
due to their special status as pseudo-Goldstone bosons. However, the pseudoscalar meson 
couplings are related to the axial current couplings, which have been analyzed in detail, 
and shown to agree with 1/A^c predictions [0. There are two couplings in the Nijmegen 
potential, the (p and a2 coupling, that deviate significantly from the 1/A^c pattern. The 
(j) meson is a pure ss state, and only couples to the nucleons through quark loops. Its 
coupling is OZI suppressed, and should be of order 1/Nc relative to the u couplings. The 
Nijmegen fit has g(j}/gLo ~ 0.1, which is a factor of three smaller than the naive 1/A^c 
prediction. The 02 coupling is even somewhat smaller. 

It must be stressed that the numerical parameters plotted in fig. 3 were obtained by 
treating the couplings as phenomenological parameters in the A^A^ potential, chosen to 
provide the best fit to A^A^ scattering data. There is no reason to assume that the A^A^ 
force is actually due to single meson exchange; in fact, the e and P contributions to the 
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Fig. 3. The couplings for the Nijmegen NN potential in ref. rescaled by fp. 
The values for this ratio predicted by large Nc QCD in eq. ( U-W o.^^ indicated by 
lines, and the shaded regions are the size of the expected 0{1/Nc) corrections to the 
leading result. The five regions in the plot (separated by vertical dashed lines) are 
(from left to right) the {I,S) = (0,0), (1,1), (1,0), (0,1), and Qfj, couplings. 



potential do not correspond to single meson exchange at all, and the 02 in the Nijmegen 
potential has a Gaussian propagator. The model subsumes such effects as 2-7r exchange, pre 
exchange, etc. within the phenomenological couplings gis- Only the pseudoscalar meson 
couplings are related to the physical meson-nucleon couplings, since the long distance part 
of the A^A^ potential is dominated by single meson exchange. Thus the agreement between 
fig. 3 and the large- A'^c prediction (|4.3|) contains more than the claim that meson-baryon 
couplings obey the It = Jt rule. We take fig. 3 to provide encouraging evidence that our 
large- A'^c analysis of the A^A^ interaction describes the qualitative features seen in nature. 



5. The Central Potential and Wigner Supermultiplet Symmetry 

It was suggested in ref. that the approximate Wigner supermultiplet symmetry 
observed in light nuclei could be explained by the 1/A^c expansion of QCD. Under the 
Wigner symmetry SU{4)w, the four nucleon states p t, p |, n t and n | transform as the 
four dimensional fundamental representation. Note that SU{4)w is distinct from the quark 
model SU{4), and that the former cannot be realized as a symmetry at the quark level. 
Nevertheless, ref. [0 argued that the 1/A^c expansion explains how SU{4)w symmetry 
could emerge as an accidental symmetry in light nuclei. As that work only examined the 
central part of the A^A^ potential, it is worth reexamining the argument. 
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Under SU{4)y\/ symmetry, a two-nucleon state transforms like 4 x 4 = 6^ + IO5, 
where the subscripts A and S denote the antisymmetric and symmetric combinations. 
Under spin x isospin, these representations decompose as 

6^ (0,1) + (1,0), 

(5.1) 

10^ (0,0) + (1,1) . 

If the two nucleons are in an even partial wave, they must be in a totally antisymmetric 
spin ® isospin state, so they are in a (0, 1) or (1, 0) state, i.e. in 6^ of SU{4)w. If the two 
nucleons are in an odd partial wave, they must be in a totally symmetric spin isospin 
state, so they are in a (0, 0) or (1, 1) state, i.e. in IO5 of SU{4)w- 

We have shown in §3 that the leading contributions to the A^A^ potential (at strength 
Nc) are Vg", and V^^. The first two terms correspond to the operators 

1, ai • 0-2 Ti • T2, (5.2) 

which have the same value on (0, 1) or (1, 0), i.e. they have the same value on the entire 
6 representation of SU{4)w- Thus at leading order in Nc, the central potential respects 
Wigner SU(4)w symmetry if the two nucleons are in an even partial wave. Violation of 
Wigner SU{4)y\/ from the central potential in the even partial waves is an 0{1/N^) effect. 
The operators eq. (|5.2| ) have different values on the (0,0) and (1,1) representations, and 
so break SU{4)w symmetry when acting on the 10 representation of SU{4)w- Thus the 
central potential breaks Wigner SU{4)w symmetry at leading order in the odd partial 
waves. The tensor force Vf also violates Wigner SU{4)w symmetry at leading order, in 
all partial waves. 

Nevertheless, there is reason to expect to see Wigner symmetry in light nuclei. The 
nucleons inside a nucleus have low momentum, so the dominant interaction is s-wave 
scattering, with higher partial waves being kinematically suppressed. Furthermore, the 
tensor mean field is small in nuclei. Therefore all of the leading order violations of SU{4)y[/ 
may be expected to be small. 

So why is SU{4)w not evident in heavy nuclei? At subleading order (a relative 
the potentials V^, Vq*^, and V^g all break the Wigner symmetry. The mean fields of 
all but the spin-orbit force are small in nuclei. However, the importance of the spin-orbit 
force grows like A^/'^, proportional to the number of particles in the maximum angular 
momentum shell. Therefore, for large A, the spin-orbit force is expected to overcome the 
suppression and destroy the approximate Wigner supermultiplet symmetry. It may 
be interesting to pursue this further, to determine at what values of A one might expect 
SU{4)w symmetry to fail. 
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6. Summary and Conclusions 

The 1/Nc expansion has been shown elsewhere to be a useful tool in analyzing the 
properties of baryons i; the analysis presented here and in ref. shows that it also provides 
a useful tool for understanding qualitative features of the nuclear force. In particular, we 
have computed the relative strengths of the various components of the A^A^ interaction in 
the 1/Nc expansion (Table 1) and argued that the predicted patterns are reproduced in 
phenomenological models of the A^A^ force (fig. 3). We also extended the argument of ref. 
that the approximate Wigner supermultiplet symmetry observed in light nuclei (see 
for examples and references) is in fact understandable in terms of the 1/A^c expansion. We 
are aware of no other explanations for this peculiar SU{4)w symmetry. 

Aside from obtaining directly from QCD a qualitative explanation for the spin, isospin 
and tensor structure of the A^A^ potential, it is hoped that the 1/N(. expansion could serve 
as a guide toward better understanding the interactions of baryons with strangeness, where 
the experimental data is much poorer. To this end we have included the three flavor analy- 
sis in Appendix A. Our hope is that this could prove useful for understanding hypernuclei, 
as well as matter in extreme conditions where strangeness may play a signiflcant role, such 



as in heavy ion collisions, or dense matter with kaon condensation [|T9| or hyperons. 
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Appendix A. Three Flavors 



The quark operator basis for three flavors is denoted by 

5^=Q^yQ, f^ = Q^T^Q, ^^^ = Q^yT^g, (A.l) 

where Q = (w, d, s) and are the creation and annihilation operators for the three 
light quark flavors, i = 1,2,3 and A = 1, ... ,8. T" are the standard SU{3) matrices 
in the fundamental representation, normalized so that trT"^T^ = 5^^ /2. These 1-quark 
operators act on a baryon state which is the completely symmetric tensor product (in spin 
® flavor) of Nc quarks. 

It is convenient to break the operator basis ( |A.1|) for the SU{6) generators by sepa- 
rating Q = (m, d, s) into q = {u, d) and s. Under this decomposition <S\ T"^ and break 
up into linear combinations of 

S'^q^-q, r^q^-q, G'-^q^—-q, 

and y*"t and -fir"^ which are the hermitian conjugates of y*" and -ftT". For baryons with 
Nc quarks, and strangeness of order one, G*" is of order N^., Y"^" and are of order 
\/iVc, and 5"^^, SI and Ns are of order one 0. Note that F*" and K"" are strangeness 
changing operators. 

The (properly normalized) SU{6) generators v^^*"^, and <S^/\^ are collec- 

tively denoted by A^. The operator basis for two and three flavors are summarized in 
Table 2. An expansion using the operator basis ( |A.2D gives us the predictions of the 1/Nc 
expansion for three flavors, without assuming SU{3) symmetry. One can also impose 
SU{3) symmetry, which places additional restrictions on the flnal result. The results for 
two flavors are obtained by using only the operators 5"*, and G*". 



TABLE 2 



^ of Flavors 




Spin 


Flavor 


Spin — Flavor 


All 


2 






fa 






3 






fA 








AS = 










3^2 


AS = 1 












AS = -1 




a 


-y>it 

a 
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The results of the paper can be generahzed to the case of three hght flavors. The 
analysis is more complicated because one also has to include operators Ng, SI, Y^"- and 
that involve the s quark. We will simply give the final results here. 

The 1/Nf, AL = interaction is 

The Ng term violates SU{3) symmetry, so its coefficient is proportional to SU{3) breaking 
in the baryon sector, which is parameterized by e, a dimensionless number of order 0.3. 
It is clear from Eq. (|A.3| ) that the Ng term violates SU (6) symmetry but respects SU (4) 
symmetry, so that SU{6) violation is of order e/Nc, but SU{4:) violation is of order 
The AL = 1 interaction for three flavors is 



^ r=0 

Nr-2 I AiArf^A I niA 



+ E ^2, 



r=0 



gr^i + GrTj \ / Ai • A 2 



r=0 



+ E ^i'-^ 

r=0 

+ E^: 



/ \ 



^ gi-^gr^L + g2^gj^>gL ^ Ai-A 2 
ATS iV2 



which is the three-flavor generalization of (|3.9| ). Time reversal and parity invariance re- 
quires the coefficients in ( p.9|) to be of the form 

Pin ^ Pout 5 

times an arbitrary function of q^, and [q- k)^ . As for the case of two flavors, the 
coefficients in eq. ( [A.4| ) are of order 1/Nc, so that the AL = 1 amplitude is of order 
relative to the central potential. 
The AL = 2 amplitude is 

where the coefficient f^^ is a symmetric traceless tensor that depends onp^^ andp^^^^. This 
is the three-flavor generalization of ( |3.11| ). 
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